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ABSTRACT 

The cold dark matter scenario predicts that a large number of dark subhalos should be located 
within the halo of each Milky-way sized galaxy. One tell-tale signature of such dark subhalos could 
be additional milliarcsecond-scale image splitting of quasars previously known to be multiply-imaged 
on arcsecond scales. Here, we estimate the image separations for the subhalo density profiles favoured 
by recent N-body simulations, and compare these to the angular resolution of both existing and 
upcoming observational facilities. We find, that the image separations produced are very sensitive 
to the exact subhalo density profile assumed, but in all cases considerably smaller than previous 
estimates based on the premise that subhalos can be approximated by singular isothermal spheres. 
Only the most optimistic subhalo models produce image separations that would be detectable with 
current technology, and many models produce image separations that will remain unresolved with 
all telescopes expected to become available in the foreseeable future. Detections of dark subhalos 
through image-splitting effects will therefore be far more challenging than currently believed, albeit 
not necessarily impossible. 

Subject headings: Dark matter - galaxies: halos - galaxies: dwarf - gravitational lensing - quasars: 
general 



1. INTRODUCTION 

While the cold dark matter (CDM) scenario has been 
very successful in explaining the f ormation of la rge-scale 
structures in the Universe (see e.g. lPrimackl2003l for a re- 
view) , its predictions on the scales of individual galaxies 
have not yet been confirmed in any convincing way. One 
particularly interesting feature of the CDM model is the 
high level of halo substructure generated. According to 
current simulations, each galaxy-mass CDM halo should 
contain large numbers of subhalos (t ypically accounting 
for k, 5-10% of its to tal mass; e.g. iNurmi et all 120061 : 
iDiemand et al.ll2007a[ ) in the dwarf-galaxy mass range. 
However, these subhalos do not appear to correspond 
to luminous structures, as the dark halo of the Milky 
Way would then contain a factor of 10-100 more satellite 
galaxies than observed, provided that each subhalo cor- 
respo nds to a luminous dwarf galaxy (e.g. iKlvpin et alJ 
1999). A similar lack of dwarf galaxies compared to the 
number of dark halos predicted is also evident on scales 
of galaxy groups (jTullv et al.ll2002t ). One way out of this 
problem is to assume that most of th ese low-mass halo s 
correspond to so-called dark galaxies (jVerde et al] f2002) , 
i.e. objects of dwarf-galaxy mass which either do not con- 
tain baryons or in which the baryons have not formed 
many stars. While a number of very faint satellite galax- 
ies have recently b een detected (e.g. IZucker et al.ll2006l ; 
iSimon et al.l [2006) , it is still far from clear that these 
exist in sufRcent n umbers to account fo r the subhalos 
predicted by CDM (jSimon fc Gehall200l . 

Gravitational lensing may in principle offer a route 
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to detecting completely dark galaxies, and the mag- 
nification associated with dark halo substructures has 
long been suspected to be the cause of the flux 
ratio anomalies observed in multiply- imaged quasar s 
(e.g. iMao fc Schneiden [19981 iKochanek fc Dalai l200l . 
However, the subhalos predicted by current CDM 
simulations appear too few to expla in the flux ra- 
tios of many of these systems (e.g. IMao et al.1 120041 : 
iMetcalfl I2005at IMaccio fc Mirandal l2006h . This in- 
dicates that other mechanisms, like l ensing by low- 
mass field halos outside the main lens (|MetcalJ [2005b: 
iMiran da & Maccio 200 7]), stellar microle nsing in the lens 
galaxy (jSchechter fc Wa mbsganss 2002) an d absorption 
or sc intillation in the interstellar medium ([Mittal et al.l 
12001 may also be at work. 

A similar but more direct way of putting the CDM 
subhalo predictions to the test could be to search for 
small-scale image separations produced b y strong lens- 
ing from subhalos. lYonehara et al.l (|2003l ) argue that by 
targeting quasars that are already known to be multiply- 
imaged on arcsecond scales, there should be a significant 
probability of detecting image splitting by subhalos at 
scales of milliarcseconds (so-called millilensing or some- 
times mesolensing). If the sources are extended and can 
be resolved at scales smaller than the Einstein radii of 
the subhalos, then not only the number densities of the 
subhalos but also their internal de nsity profiles can be 
constrained (jlnoue fc Chiba|[2005bl ). It has been argued 
that this techniqu e may become obse rvationally feasi- 
ble with ALMA (jlnoue fc C hiba 20 05al) or fu ture space- 
VLBI missions like VSOP-2 (|Inoue fc CInball2005cl ). If 
correct, this would constitute a major step forward in 
the study of dark halo substructures. 

These proposals do however rely on the assumption 
that the image separations (or Einstein radii) of subhalo 
lenses are similar to those produced by singular isother- 
mal spheres. This is an assumption that is difficult to 
justify, since neither simulations, theoretical arguments 
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nor observations favour a density profile of this form for 
dark matter halos in the relevant mass range. 

N-body simulations based on the CDM paradigm typi- 
cally predict dark matter halos to have inner density pro- 
files of the fo rm p(r) oc r~ a with central den sity slopes 
a « 1 (e.g. iNavarro. Frenk. fe White! [1996, hereafter 
NFW), whereas the singular isothermal sphere (SIS) as- 
sumes a — 2. While the SIS has proved to be a successful 
model for observed galaxy-m ass lenses (e.g. iRusin et al.l 
120031 : iKoopmans el al.ll2006f ). this is believed to be due 
to the luminous baryons residing in the inner regions of 
these objects. This component contributes substantially 
to the overall mass density in the centre, and its forma- 
tion over cosmological time scales may also have caused 
the CDM halo itself to contract, the reby steepening the 
inner slope of its density profile (e.g. jGnedin et alj|2004t 
Maccio et al]|2006t iGustafsson et all2006l:lKampakogloul 



2006). 



Dwarf galaxies are well-known to be more dark mat- 
ter dominated than their more massive counterparts, and 
some of the mechanisms proposed for keeping the miss- 
ing satellites of the Local Group dark assume that most 
of their baryons have been lost early in the history of 
the Universe, or are prevented from collapsing into their 
centr al regions (e.g. iBarkana fe Loeblfl999t iRead et all 
120061 ). Hence, it seems reasonable that the inner slope 
a of their density profi l es sho uld be shallower than that 
of an SIS. Indeed, [Mi ([20031 ) finds that to explain the 
distribution of angular separations among strong lenses, 
the lens population must shift from an SIS-type profile 
(a = 2) to something more resembling the NFW profile 
(a w 1) at halo masses of < 10 12 M©. Observational con- 
straints on the slope of the innermost density profile of 
both field and satellite galaxies in the dwarf galaxy mass- 
range typically favour a slope even sha l lower than this 
(a » 0-1; e.g. Ide Blok fc Bosmal [2002t ISpgkkens et al.l 
mM iZackrisson et al.l [2006: Gi lmore et all 12007ft . This 
is potentially bad news for the detectability of low-mass 
subhalos through millilcnsing effects, since the image sep- 
aration expected from strong lensing by a halo of fixed 
virial mass is extremely sensitive to the slope of its inner 
densit y profile. In fact, as pointed out bv lTrentham et all 
(|2001h . the image separation expected from strong lens- 
ing by dark galaxies would be unobservably small if their 
density profiles are of NFW type. An inner slope of the 
density profile even shallower than this would just make 
the problem worse. 

As low-mass halos in the field fall into the potential 
wells of larger halos, tidal interactions with their host 
halos and with other subhalos will cause these objects 
to lose mass. Since this mass loss preferentially takes 
place in the outer regions of the subhalos, their den- 
sity profiles change ove r time (e.g. lHavashi et all 120031 : 
iKazantzidis et aiT [2004). thereby altering their lensing 



properties compared to more isolated field halos of the 
same mass. Here, we take a critical look at the prospects 
for strong-lensing detections of dark subhalos in the 
dwarf-galaxy mass range, by deriving the image sepa- 
rations expected for subhalo density profiles favoured by 
current simulations. 

2. DARK HALO DENSITY PROFILES 

When deriving the subhalo density profiles and the 
corresponding image separations, we assume a ACDM 



cosmology with fl\ = 0.7, Qm — 0.3 and h = 0.7 
(Hq = lOO/i km s _1 Mpc -1 ). As we are mainly inter- 
ested in image-splitting of the individual macroimages of 
multiply- imaged quasars, we compute all density profiles 
at a typical lens redshift of z\ = 0.5 and adopt a source 
redshift of z s = 2.0. The impact of of other choices for 
z\ and z s on our results are discussed in section 3. We 
moreover assume that halos contain no material beyond 
their virial radius r v i r (zi), defined as the radius at which 
the mean enclosed density equals A V j r times the criti- 
cal density of the Universe at redshif t z\. Here, A v j r is 
calcul ated using approximations from iBrvan fc Normanl 
(fl998h . 

The majority of papers on the density profiles of CDM 
halos have so far focused on relatively isolated halos. 
When such objects are accreted by more massive halos 
and become subhalos, substantial mass loss occurs, pref- 
erentially from their outer regions. Because of this, some 
sort of modification of these density profiles is in order 
when applying them to the study of subhalos. However, 
since some authors have previously modelled the lensing 
properties of subhalos using the original, unstripped pro- 
files, we will here consider such models as well, so that the 
the importance of the stripping on the resulting image 
separations becomes clear. It is moreover possible that 
the original profiles may be reasonable approximations 
for recently accreted objects or low-mass halos which are 
spatially correlated with the parent halo but formally 
located outside its vir ial radius (e.g. iNurmi et al.1 120061 
iDiemand et al.ll2007bh . Even though the latter objects 
are located far away from the central regions of their host 
halos, they are predicted to exist in great numbers and 
may end up projected close to the images of strongly- 
lensed quasars. 

2.1. Unstripped profiles 

The original NFW density profile, relevant for rela- 
tively isolated CDM halos, is given by: 

Pi 



(r-/r s ,NFw)(l + r/r s ^ FW ) 2 ' 



(1) 



where rs,NFW is the characteristic scale radius of the 
halo. 

Because of the persisting controversies concerning the 
most appropriate fitting formula fo r the density profiles 
of CDM halos, we also consider the lMoore et al.l (|1999bl 
hereafter M99) profile, which - while having an inner 
density profile steeper than that of the original NFW 
formula - gives a fit of c omparable quality to simulation 
data for low-mass halos ([Navarro et al.l l2004): 



po 



(2) 



(r/r SM 99y- 5 [1 + (rAs.Mgg) 1 ' 5 ] ' 

INavarro et alj (|2004l hereafter N04) themselves do how- 
ever advocate yet another profile, 

PmmM = ps,N04 exp(-2//3 [(r/rs^o^ 13 - l] ), (3) 

with the property of asymptotically approaching a finite 
central density. Here, ps,noa is the density at rs,N04j i-c. 
Pi/4, and (3 « 0.17 for halos in the dwarf-galaxy mass 
range. 

To derive the R s parameters of the NFW, M99 and 
N04 profiles, we compute these in units of the radius 
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r_2 at which the power-law slope of the density profile 
obeys dlnp/dlnr = —2. This quantity is linked to the 
virial radius r V j r through the conc entration paramete r 
c = f* V i r /r_2 7 with a median value of (Bull ock et a.ll206"lf) : 



1 



Mh 



1.5 x 10 13 M„ 



-0.13 



(4) 



Contrary to many other c(M, z) models, this one has 
been d emonstrated to ex tend into the dwarf-galaxy mass 
range (jColm et al.ll2004 L Once the rs parameters are de- 
termined, pi, ps and po can then be derived by requiring 
that the correct mass is contained within r v i r . 

In contrast to these CDM-motivated halo models, the 
commonly used SIS model has a density profile of the 
form: 

PS1S = 2^' (5) 
where o~ v is the line-of-sight velocity dispersion. 

2.2. Stripped profiles 

Based on N-body simulations of subhalos orbit- 
ing within a stat ic and spherical host potential, 
lHavashi et alj (|2003l hereafter H03) suggested a subhalo 
density profile that represents a simple modification of 
the NFW profile of its progenitor: 



ft 



PH03(r) = -7—1—1 ^PNFWM 

1 + {r/r tc y 



(6) 



where ft, the central density reduction factor, and r te , 
the effective tidal radius of the subhalo (in units of 
the progenitor scale radius rs), are both related to the 
mass fraction of the subhalo that remains bound, mt>nd, 
through relations given in H03. For simplicity, we calcu- 
late the properties of the subhalo progenitor profile by 
assigning it an approximate mass M prog = M su b/mbnd- 

Based on similar sim ulations at higher resolution, 
iKazantzidis et ail (|2004l hereafter K04) et al. instead 
find a density profile of the form: 



PK0i(r) = kr 7 exp( ), 



(7) 



where 7 w 1 ± 0.05, 0.1 < r b /r S NFW < 1.2 and rg !NFW 
refers to the corresponding parameter of the progenitor 
halo, which we assign a mass analogous to the H03 case. 
Here, k is a normalization factor which we introduce to 
scale the subhalo to its current mass mbnd-Mprog- In the 
following, we adopt 7 = 1 and rb/rs,NFW = 0.75. 

In the case of both the H03 and K04 profile , we assume 
the subhalo mass to be entirely contained within r v ; r of 
its progenitor and adopt rn b nd = 0. 03 as a typical value 
for su bhalos of galaxy-mass halos (|van den Bosch et al.l 
I2005T) . 

While many previous investigations of the lensing ef- 
fects produced by subhalos have assumed a sharp den- 
sity cutoff at some outer radius, it is clear from the H03 
and K04 profiles that the density decrease in the outer 
regions is expected to be gradual. To demonstrate the 
effects of this, we consider two different modifications of 
the M99 and N04 profiles - one with a sharp truncation 
radius r t such that the mass contained inside becomes 
M su h(< ft) = m-bnd-Mprog, and one in which the halo is 



gradually stripped of material according to the function 
proposed by H03: 

ft 

PM99/N04,strippod( r ) = Y+Jrjr ) 3 ^ M99 / N04 ( r )" ^ 

Lacking a better prescription, we adopt the /t(?7ibnd) and 
rte(wbnd) relations presented in H03, although these will 
in practice depend on the details of the progenitor halo 
model assumed. 

When considering the impact of stripping and/or trun- 
cation on lensing properties, we have chosen to compare 
halos of fixed mass, which is defined as the virial mass 
for unstripped halos, and as the total subhalo mass after 
stripping for the subhalos. This choice - which reflects 
the tradition in the field of subhalo simulations, where 
one usually present the mass function of the subhalos 
themselves and not their progenitors - has important 
consequences. For instance, a sharply truncated halo of 
a certain mass has a much smaller spatial extent than an 
untruncated halo of the same mass. Consequently, the 
density in the inner regions of the truncated halo is nec- 
essarily higher than that in the untruncated halo of the 
same mass, and therefore, the truncated halo will have a 
larger Einstein radius. However, a more gradual density 
decrease in the outskirts can in principle drive the Ein- 
stein radius in either direction, depending on the details 
of the stripping process. As it turns out, the stripping 
procedure adopted here (equations ^ and ©) tends to 
produce Einstein radii that are smaller than those of the 
unstripped profiles. 

3. STRONG-LENSING IMAGE SEPARATION VERSUS 
OBSERVATIONAL ANGULAR RESOLUTION 

To first order, the image separation produced by an ex- 
tended object with a density that decreases as a function 
of distance from the centre is given by 



A0 » 2R, E /D oh 



(9) 



where D a \ represents the angular-size distance between 
observer and lens, and Re represents the linear Einstein 
radius. The latter is defined as_the radius inside which 
the mean surface mass density E of the lens equals the 
critical surface mass density 



E(< R E ) = E c = 



c 2 D n 



4ttGD q iDi s ' 



(10) 



where D os and D\ s are the angular-size distances between 
observer and source, and lens and source, respectively. 

From the density profiles described in the previous sec- 
tion, we numerically compute the mean surface mass den- 
sity profiles E(< R), shown in Fig. QJi and b in the case 
of a subhalo lens mass of 10 10 M Q . At the adopted lens 
rcdshift, the smallest linear radius plotted corresponds 
to an angular Einstein radius of w 3 x 10 -6 arcseconds, 
which is below the angular resolution of any telescope 
that will become operational in the foreseeable future. 

In fact, only the SIS and M99 profiles cross the E c di- 
vide (gray horizontal line) at radii substantially larger 
than this, which means that only these profiles will give 
rise to detectable image separations for a subhalo_of this 
mass. In its original form, the NFW profile gives E(< R) 
lower than E c _at these radii, whereas the truncated ver- 
sion leads to E(< R) above the threshold. Both the H03 
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Fig. 1. — a) Mean surface mass density £(< r) profiles for a 10 10 Mq subhalo at z\ = 0.5: SIS (thin solid), truncated SIS (thick solid), 
NFW (thin dash-dotted), truncated NFW (thick dash-dotted), H03 (thick dotted) and K04 (thin dotted), b) Same for N04 (solid with 
triangles), stripped N04 (solid with squares), truncated N04 (solid with cirles), M99 (medium dashed), stripped M99 (thin dashed) and 
truncated M99 (thick dashed) profiles. In both panels, the horizontal gray line indicates the critical surface mass density for strong lensing, 
S c . The truncated models have been constructed from density profiles for halos of mass 10 10 mr\ Mq (with m^j = 0.03), abruptly 
terminated at a radius enclosing a mass of 10 10 Mq . The stripped versions of the N04 and M99 profiles have been constructed using JHJ. 



and K04 subhalo models, which should be far more realis- 
tic than any sharp truncation, do however predict 
smaller than that of an undisturbed NFW halo of the 
same mass. The N04 profile (in both its original, stripped 
and truncated form) also ends up below the threshold, 
due to its finite-density core. The effects of truncation 
versus gradual stripping are qualitatively similar for the 
NFW, N04 and M99 halos - a sharp truncation increases 
the Einstein radius whereas a gradual stripping decreases 
it relative to the Einstein radius produced by the original 
profile. This indicates that previous investigations based 
on sharp outer truncations are likely to have substantially 
overestimated the image separations. 

In Fig. [21 we plot the image separations predicted for 
10 4 -10 11 Mq subhalos against the angular resolution 
of a number of planned or existing observational facil- 
ities, operating at a wide range of wavelengths. These 
include the proposed MAXIM pathfind er in X-rays 6 ; 
VLTI with the proposed VSI instrument (Malbct ct al. 
I2006D . the planned GALA 7 and SIM PlanetQuest 8 satel- 
lites in the optical/near-infrared; the currently avail- 
able EVN 9 , HSA 10 , VLB A 11 arrays plus the planned 
ALMA 12 , EVLA 13 , SKA 14 arrays, and also space- VLBI 
with the planned VSOP-2 15 satellite at radio wave- 
lengths. Please note that here we consider only the 
best resolution limits attainable with these telescopes, 
whereas the resolution at the wavelengths that maximize 
the number of observable high-redshift sources may be 
considerably worse. 

It is immediately evident from Fig. [2] that there are 
large differences between the image separation predic- 

6 maxim. gsfc . nasa.gov / pathfinder . html 

7 www.rssd.esa.int/index.php?project=GAIA 

8 planetquest.jpl.nasa.gov/SIM/ 

9 www.evlbi.org 

10 www.nrao.edu/HSA/ 

11 www.vlba.nrao.edu 

12 www.alma.info 

13 www.aoc.nrao.edu/evla/ 

14 www.skatelescope.org 

15 www.vsop.isas.ac.jp/vsop2/ 



tions of the various halo models. As the discrepancy 
between the number densities of luminous galaxies and 
dark matter halos does not start to become severe until 
the halo mass drops below < 10 10 Mq (e.g. IVerde et al.l 
2002; Ivan den Bosch et aTll2003h . subhalos at masses be- 
low this limit need to produce measurable image sepa- 
rations (fl>4x arcsec for VSOP-2, which has the 
best theoretical resolution among the telescopes included 
in Fig. [2]) in order for dark galaxies to be detectable 
through image-splitting effects. Out of the halo models 
tested, only two actually meet this criterion without ad- 
hering to sharp truncations: the SIS and the M99 halos. 
The H03 and K04 profiles both give image separations 
smaller than 10 -6 arcsec for all the halo masses con- 
sidered and are therefore completely outside the plotted 
region. Even in the optimistic case of an M99 halo (in 
either its original or stripped form, whereas the sharp 
truncation, as discussed previously, is not considered re- 
alistic), the image separations are a factor of w 3-7 
smaller than those predicted for a SIS (and « 30-60 
times smaller than those of a truncated SIS), render- 
ing only the few most massive subhalos (~ 10 10 M Q or 
slightly higher) detectable at ~ 0.01" resolution (GAIA, 
SIM and ALMA). At milliarcsecond resolution (VLTI 
and SKA), dark galaxies with masses > 10 9 M© may 
become detectable. To probe further down the sub- 
halo mass function, submilliarcsecond-resolution facili- 
ties (MAXIM pathfinder, HSA, EVN, VLBA or VSOP-2) 
will be required. 

These image separations have been computed for fidu- 
cial lens and source redshifts of z\ = 0.5 and z s = 2.0, but 
the overall picture does not change substantially for other 
realistic choices of these parameters. A higher z s implies 
a lower critical surface mass density S c for certain z\, 
which for a fixed subhalo density profile will increase the 
Einstein radius and boost the image separation A9. This 
boost becomes more pronounced for objects with flat in- 
ner E(< R) profiles, like the NFW, N04, H03 and K04 
density profiles (see Fig. [J). While the density profiles 
themselves also change with z\ through the redshift evo- 
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Fig. 2. — Subhalo mass versus image separation A8 for those density profiles from Fig[T]that give rise to image separations on scales larger 
than microarcseconds. The angular resolution of a number of existing and planned observational facilities have been indicated by horizontal 
dashed lines, marked with labels (see main text for details), a) The different diagonal lines represent SIS (thin solid), truncated SIS (thick 
solid), M99 (medium dashed), stripped M99 (thin dashed) and truncated M99 (thick dashed) subhalo models, b) The different diagonal 
lines represent NFW (thin dash-dotted), truncated NFW (thick dash-dotted) and truncated N04 (solid with circles) subhalo models. 



lution of the virial radius and the concentration param- 
eter, this has a minor effect on AO, compared to impact 
of the S c (zi,z s ) dependence. We find, that by allowing 
source redshifts as high as z s = 5.0, the image separation 
of a 10 10 Mq NFW-type subhalo can be increased by a 
factor of up to « 30. Even though this boost factor may 
seem large, the resulting image separation is still only 
AO ~ 10 -6 arcseconds, and hence far below the observ- 
able range. For an M99 profile, the corresponding boost 
is only w 30%. Hence, subhalos need to have steep in- 
ner density profiles (of approximately the M99 type) to 
become observable, even if a sample of multiply- imaged 
quasars with lens and source redshifts chosen to maxi- 
mize the subhalo image separations would become avail- 
able. 

4. EXTERNAL CONVERGENCE AND SHEAR 

The AO estimates presented in the previous section are 
based on the assumption that the subhalos can be treated 
as isolated objects. The effects of external convergence 
K and shear 7 due to the galaxy (and hal o) hosting these 
subha los have thereby been neglected. lYonehara et al.l 
( 2003) explored the consequeces of non-zero k and 7 in 
the case of SIS subhalos and found the external poten- 
tial to have a non-negligible impact on the strong lensing 
properties of such objects. However, their study sug- 
gests that the resulting image separations A0 Kn are typ- 
ically within a factor of w 3 from those derived in the 
case of k = 7 = 0. If applied to the predictions plot- 
ted in Fig. [2j a boost factor of this magnitude would 
be insufficient to challenge our main results, namely 
that only subhalos with inner density profiles as steep 
as those of M99 (or SIS) models produce image sep- 
arations that can be resolved with current or planned 
telescopes. The question remains whether the boost fac- 
tor /boost = A6' Kj7 /A6' k= o,7=o can become considerably 
larger for some of the more realistic density profiles con- 
sidered here. 

To investigate this, we use ray-tracing simulations to 
numerically assess the distribution of /boost- For the 



macrolens (i.e. the halo and galaxy hosting the subhalo) 
we adopt an SIS density profile with a v = 150 km s _1 . 
At z\ = 0.5 and z s = 2.0, the corresponds to a linear (an- 
gular) Einstein radius of 2.5 kpc (0.4 arcsec). Subhalos 
are then distributed within this structure, assuming an 
NFW profile with c = 10 for the subhalo component of 
the CDM in the macrolens. Sinc e the optical depth fo r 
image splitting by subhalos is low ijYonehara et al.l l2003) , 
the /boost distribution is not very sensitive to the exact 
spatial distribution of the subhalos within the host halo. 
For simplicity, we assume all subhalos to have the same 
mass in each simulation, but repeat the simulations for 
different subhalo masses to explore the dependence of 
/boost on the subhalo mass. The sources are assumed 
to be point-like and are distributed on a regular grid in 
the source plane. To assess the effects of magnification 
bias, two different cases are considered: one in which 
no magnification threshold is imposed, and one in which 
only source positions which produce total magnifications 
fi > 10 are analyzed. Source positions for which multiple 
images are not produced are always rejected. 

The vast majority of the resulting macroimages turn 
out to be unaffected by the subhalos due to the low op- 
tical depth, and are therefore discarded. The proper- 
ties of the Einstein ring of the subhalo is calculated for 
the remaining macroimages, and the image separation 
A# Kj7 estimated from its angular diameter. In the case 
of 7 =/= 0, the Einstein ring becomes an ellipse, and A0 Kn 
is then estimated along the major axis. This leads to a 
systematic overestimate of A# K)7 , which ensures that the 
resulting boost factors are conservative upper limits. 

In the case of no magnification threshold, the resulting 
distribution of /boost is shown in Fig. [3] for an M99 sub- 
halo (solid line) and an H03 subhalo (dashed) , both with 
Mmb = 10 10 Mq. As discussed in Section 3, these sub- 
halo models span the range from effectively undetectable 
(H03) to favourable for detection (M99). In both cases, 
the tail of the /boost distribution extends up to very high 
values, which means that it is in principle possible to 
find macroimages for which subhalos of masses much 
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Fig. 3. — The cumulative probability distribution P(> /boost) 
of image separation boost factors /boost due to the external poten- 
tial in the case of M99 (solid) and H03 (dashed) subhalo density 
profiles. In both cases, the mass of the subhalo is 10 10 Mq. 



lower than indicated by Fig. [2] can be detected through 
image-splitting effects. However, such macroimages are 
exceedingly rare. The probability of having /boost > 100 
is ~ 10~ 3 in the case of M99 and rj 0.03 in the case 
of H03. The expectation value for the boost factor is 
(/boost) « 2.3 for the M99 subhalo and (/boost) ~ 14 
for the H03 subhalo. This is enough to shift the AO 
prediction for the M99 subhalos fairly close to that of 
an untruncated SIS in Fig. [21 but insufficient to move 
the H03 prediction (Ar? K= o,7=o ~ 10~ 15 arcsec at this 
mass) into the detectable range. In the case of M99 
subhalos, we have also explored the mass dependence of 
(/boost), finding (/boost) « 3.1 for Af sub = 10 8 M Q and 
(/boost) w 3.6 for M sub = 10 6 M . This trend with in- 
creasing (/boost) for decreasing subhalo mass translates 
into a slight change in slope for the A#(M BU b) predic- 
tion in Fig. [21 but has no dramatic effects on our overall 
conclusions. 

Flux-limited observations are affected by a magnifi- 
cation bias, since flux threshold of the survey tends to 
favour highly magnified quasars from the vast popula- 
tion of intrinsically faint objects over modestly magnified 
(or even demagnified) quasars from the relatively small 
population of intrinsically bright objects. The exact im- 
pact of this effect depends on the flux threshold and the 
quasar luminosity function (which, in turn, depends on 
the wavelength at which the quasars are observed), and 
is outside the scope of this paper. Here, we instead make 
a crude assessment of the likely impact of magnification 
bias by rejecting all systems which do not fulfill the cri- 
terion fj, > 10. As expected, this increases the boost 
of the image separation produced by the subhalos. In 
the case of M su b = 10 10 M Q subhalos, the boost reaches 
(/boost) ~ 6.7 for an M99 profile and (/boost) « 85 for an 
H03 profile, i.e. an increase by a factor of w 3 (M99) and 
w 6 (H03) compared to the case without magnification 
bias. While this brightens the prospects of detecting sub- 
halos with M99 profiles, H03 subhalos still produce image 
separations that are too small to be resolved with any of 
the telescopes considered in Fig. [21 

5. DISCUSSION 



Our results indicate that detection of subhalos through 
gravitational image-splitting is likely to be considerably 
more challenging than suggested in previous studies, due 
to the smaller image separations predicted for subhalo 
density profiles more realistic than the SIS models often 
adopted. In fact, no currently planned telescope will 
be able to resolve the image separations produced by 
subhalos with density profiles of the type suggested by 
the most realistic simulations currently available (H03 & 
K04). We stress, however, that these simulations do not 
necessarily represent the final word on this issue. NFW 
density profiles were adopted for the subhalo progenitors 
in both the H03 and K04 simulations, which moans that 
the slope of the central density profile prior to stripping 
was assumed rather than derived from the simulations 
themselves. If the subhalos would have steeper central 
density slopes (e.g. of M99 type), these would give rise 
to image separations that could be resolved even with 
existing telescopes. 

Despite the somewhat bleak detection prospects pre- 
sented here, there are at least two effects that can po- 
tentially improve the detectability of image splitting by 
subhalos: baryon cooling and the presence of intermedi- 
ate mass black holes. 

Baryon cooling would cause t he subhalo (or its 
progenitor halo) to contract (e.g. [Gncd n c t al.l 120041; 
Maccio et a l. 2006; Gustafss on et al.f2006l : lKam Dakogrou 



2006), thereby increasing the central density and boost 
ing the image separation. It is, however, not clear how 
strong this effect is likely to be, since this depends on 
the details of the mechanism that prevents the dark sub- 
halo from forming stars. Baryon cooling is usually as- 
sumed to be associated with star formation, and most 
attempts to explain the lack of bright subhalos there- 
fore propose that the baryons have either been lost early 
in the history of the Universe, or are prevented from 
cooling by th e ultraviolet backgroun d provided by reion - 
ization (e.g. iBarkanafc Loebl 119991 : iRead et al.l l2006f) . 
Neither mechanism is likely to result in any significant 
halo contraction. Nonetheless, claims of dark baryons 
in the form of cold gas in galactic disks have been made 
()Bournaud et al.ll2007f ). implying that there may be some 
route for gas to cool without form ing stars or being de- 
tected by the usual H2 trac ers (e.g. lPfenniger et al"1ll994t 
iPfenniger k Combesll994r) . One may therefore speculate 
that there could be alternative solutions to the missing- 
satellite problem, in which subhalos are kept dark even 
though baryon cooling has taken place. 

Intermediate mas s black holes (IMBHs, see e.g. 
Ivan der Marel l200i iZhao fc Silkl 120051 iNovola et all 
120081 ) with masses of ~ 10 2 -10 4 M would also boost 
the image separations, if present in the centres of dark 
subhalos. This could for instance be the case if the empir- 
ical relations between the ma ss of a superma ssive black 
hole and its dark matter halo (|Ferrare se 2002) would ex- 
tend into the dwarf-galaxy mass range. Even in the case 
of k = 7 = 0, an IMBH of mass 10 4 Mq, would give 
an image separation of ps 4 x 10 -4 arcsec for z\ = 0.5 
and z s = 2.0. This is sufficient to allow VSOP-2 to re- 
solve the image splitting, regardless of the density profile 
of the subhalo hosting the IMBH. Of course, even if the 
Ferrarese relations would hold for luminous dwarf galax- 
ies, the do not necessarily do so for dark ones, since this 
depends on the formation details of IMBHs. Moreover, 
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a population of halo IMBHs not associated with subha- 
los could form an undcsired background of millilensing 
events that would obfuscate attempts to study subhalos 
through image-splitting effects. 

Finally, there is an important issue related to the limi- 
tation of current CDM halo simulations. Current N-body 
simulations only resolve scales down to R ~ 0.001r v i r , 
whereas for the steeper density profiles - i.e. those 
producing detectable image separations for dwarf-galaxy 
masses - a non- negligible fraction of £(< r) at the Ein- 
stein radius comes from radii smaller than this. This 
makes the exercise of deriving dark halo lensing prop- 
erties from fitting functions based on current N-body 
data somewhat hazardous. Unless theoretical arguments 
(e.g. lHansen fc Stad el 2006) can be used to determine the 
slope of the density profile at even smaller radii, robust 
predictions for the image separations produced by low- 
mass halos may be very difficult to derive. Observational 
limits on the lensing optical depth as a function of A8, 



based on searches for subhalo image splitting, may even- 
tually turn out to be a useful way of setting constraints 
on the internal structure of dark matter halos at scales 
below the resolution of numerical simulations. Of course, 
the optical depth does not depend only on the density 
profile of the subhalos, but also on their mass function 
and the spatial distribution within the macrolens. De- 
tailed estimates of the optical depth for strong lensing by 
subhalos in the dwarf-galaxy mass range is outside the 
scope of this study, but will be covered in paper 11 of this 
series (Riehm et al. 2008, in preparation). 
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